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The large-angle convergent-beam electron-diffraction technique, available in a conventional transmission 
electron microscope, has been employed in order to characterize the defective layer present in heavily damaged 
high-dose hydrogen-implanted and low-temperature annealed (7^500 °C) single-crystal silicon in terms of 
lattice strain and static disorder. These quantities have been measured to determine the mean relaxation volume 
and atomic concentration of defect clusters present in the implanted layers in order to investigate the structural 
features causing the reverse annealing phenomena observed in ion channeling measurements. In particular, the 
mean relaxation volume detected in the 300 °C 2-h annealed sample (0.1 nm 3 ) results in a factor three times 
higher than that measured in the as-implanted sample; on the contrary the mean atomic concentration of 
clusters does not vary appreciably after this thermal treatment. This experimental evidence suggests a noncon- 
servative growth of clusters in the low-temperature annealing regime. After annealing at 500 °C for 2h, an 
increase of the relaxation volume and a significant decrease of the mean concentration is found, thus suggesting 
that only after this thermal treatment, producing in the meantime intrinsic defects with extended internal 
surfaces, defects seem to follow a conservative ripening. 
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I. INTRODUCTION 

Hydrogen in silicon has been the subject of intense study 
as it is known to produce defect passivation, deactivate dop- 
ants, and induce point or extended defects. 1 In high concen- 
trations, introduced by in-diffusion l or by ion 
implantation, 2,3 hydrogen is known to produce extended pla- 
nar defects (platelets), having {111} and, mainly, {100} habit 
planes, evolving in large "lattice cracks" observed in 
samples implanted at ion doses of the order of 10 17 cm" 2 
and annealed at r>400°C. Platelets are then believed to 
trap H 2 molecules in a stressed state, building internal pres- 
sure that ultimately leads to silicon cracking. This phenom- 
enon is the basis of a silicon-on-insulator technology known 
as a "smart cut." 4 

The damage evolution in samples implanted at doses up 
to 2X10 16 cnT 2 has been studied using Rutherford back- 
scattering spectrometry and channeling (RBSC), elastic re- 
coil detection analysis, and transmission electron microscopy 
(TEM). 2 It has been found that most of the interstitial-type 
damage created by the ion implantation process anneals out 
after 2 h of thermal treatment at 7=200 °C. Samples iso- 
chronal annealed at 200<7<400°C for 2 h reach a quasis- 
tationary state characterized by a negligible amount of hy- 
drogen loss. In these samples, at a depth corresponding to the 
hydrogen projected range, the raise of a displacement field is 
detected in RBSC spectra and platelets are observed in TEM 
images. Hydrogen effusion starts for 400<7<500°C, in 
this temperature range, the dimension and concentration of 
extended defects remain almost constant whereas RBSC 
shows an dechanneling maximum at r=400°C. The com- 
parison between RBSC and TEM measurements have 
strongly suggested that extended defects detected by TEM 
were not the primary source of the observed dechanneling in 



the low-temperature annealed samples. 2,5 The characteristic 
"reverse annealing" observed in RBSC measurements has 
been ascribed to small clusters of H-vacancy complexes 
hosting H 2 molecules, too small to be observable by TEM. 
At present there is no structural characterization of these 
small defects, probably because of the experimental difficul- 
ties connected to the strong state of disorder present in the 
samples implanted at a high (> 1 X 10 16 cm" 2 ) hydrogen 
dose. 6 

Fourier transform infrared (IR) spectroscopy experiments 7 
have supported the idea of hydrogen trapping in intrinsic 
defects. This study has shown that isochronal annealing for T 
between 200 °C and 400 °C, where the amount of effused 
hydrogen is negligible, transforms the IR-active silanic hy- 
drogen bonded in multivacancy defects, present in the as 
implanted sample, in IR-active silanic hydrogen in single 
vacancy complexes (VH 3 and VH 4 ), and IR-mute hydrogen 
(likely H 2 ). Further annealing till 7>400°C allows to ob- 
serve IR peaks ascribed to {100} and {111} internal hydrogen 
terminated surfaces. 

Recent theoretical studies 8 have supported all the previ- 
ously reported experimental findings demonstrating that va- 
cancies play a central role in the nucleation and growth of 
hydrogen-related extended defects. In particular, planar ar- 
rays of second-neighbor VH 4 's seem to be the seed for both 
small bubbles and platelet growth, which is achieved through 
the formation of damage clusters involving the absorption of 
diffusing molecular hydrogen and emission of silicon inter- 
stitial atoms. 

The structural analysis of the defective layer formed after 
annealing of high-dose hydrogen-implanted silicon samples 
seems therefore, crucial for a deep characterization of the 
Si(H) system. 

The aim of the present paper is the study of the structural 
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features associated with small defects present in the defective 
Si(H) layers, responsible for the observed dechanneling be- 
havior in samples annealed at temperature. 7*^500 °C. In 
particular, the relaxation volume and the atomic concentra- 
tion of small clusters of damage will be determined. This 
goal will be pursued using large-angle convergent-beam 
electron-diffraction (LACBED) measurements of lattice 
strain and static disorder, as it is known that the atomic dis- 
placements produced by small clusters affect both the angu- 
lar position (strain) and intensity (static disorder) of dif- 
fracted beams. 9,10 

n. SAMPLE PREPARATION AND ANALYSIS 

Samples have been prepared by ion implantation of H 2 + 
into single-crystal silicon. Silicon Czochralski grown disks 
were 6 in. in diameter, (100) oriented, and /?-type doped with 
a resistivity between 30-50 ft cm. During implantation, the 
sample was tiled by 7° with respect to the beam. Wafers have 
been implanted with H 2 + at an energy of 31 keV with a dose 
of 8 X 10*" 15 cm -2 , but because of the H 2 + fragmentation at 
the surface and due to the low density of collisions, this 
implantation is equivalent to the implantation of atomic H at 
a double fluence, L6X 10 16 cm" 2 , and at an energy of 15.5 
keV. During the implantation process the samples were kept 
at liquid-nitrogen temperature; after implantation some 
specimens underwent vacuum annealing at 300 °C and 
500 °C for 2 h in order to study the thermal evolution of the 
crystal damage. 

[1 TO] TEM cross sections were prepared by mechanical 
dimpling, for a residual thickness of about 10 fjm, and low 
angle ion milling with Ar + ions at 4keV. TEM images were 
taken with a Philips EM400 operating at 100 keV. LACBED 
patterns have been digitally recorded in the nanoprobe mode 
(nominal spot size of 10 nm) by means of a 1-in. 694 Gatan 
charge couple device slow scan camera mounted on the same 
microscope. These patterns are obtained using a probe with a 
large convergence angle, which is focused in a plane slightly 
below or above the sample, thus forming a diffraction pattern 
(instead of an intermediate image) in the plane of the se- 
lected area aperture. This aperture is used to isolate one re- 
flection to prevent overlapping of the diffracted beams in the 
final detection plane. 11 Bright field LACBED (BFLACBED) 
or dark field LACBED (DFLACBED) patterns can then be 
obtained by selecting the transmitted beam or a diffracted 
beam, respectively (Fig. 1). In these patterns different re- 
gions of the sample contribute to different parts of the same 
diffraction pattern, thus allowing ah accurate comparison be- 
tween the diffracted intensity coming from different areas of 
the sample with a spatial resolution of the order of the probe 
size. This peculiarity of the LACBED technique applied to 
TEM cross sections allowed to determine the static disorder 
depth profile in H, as implanted Si single-crystal samples 12 
and both the size-induced tetragonal strain 13 and static 
disorder 14 in SiGe/Si heterostructures. 

Static disorder measurements on TEM cross sections have 
been performed by tilting the sample of about 10° from the 
[HO] zone axis around the [110] direction till the [441] 
zone axis is reached, which is a "sparse" one with sheets of 




FIG. I. Experimental LACBED geometry used in static disorder 
measurements. An electron beam with a very high convergence im- 
pinges on the sample oriented along the [441] zone axis and is 
focused below the specimen plane. The selected area aperture lo- 
cated in the corresponding plane of the microscope allow to select 
the (880) diffracted beam. The (880) dark field LACBED 
(DFLACBED) pattern is formed in the recording plane. 

atomic {hhO} planes parallel to the zone direction. Then the 
£ = 880 diffracted beam has been brought into the Bragg po- 
sition in order to' be analyzed (Fig. 1). Due to the g-r=0 
invisibility rule of the kinematical diffraction theory, the un- 
desired effect on diffracted intensities due to sample prepa- 
ration can be reasonably minimized, as (hhO) reflections are 
perpendicular to the [ 1 F0] direction, normal to the free sur- 
face of the TEM cross section. 

Strain measurements, starting again from the [1 F0] zone 
axis, have been performed on BFLACBED patterns obtained 
by tilting the specimen to about 3.5° around the [110] direc- 
tion until the Bragg contour of {333} crystallographic planes 
crosses the defective region, avoiding intersection with the 
Bragg contours corresponding to different reflections (see 
Fig. 2). 

m. LACBED STATIC DISORDER AND STRAIN 
DETERMINATION 

A. Static disorder measurements 

The information on static disorder is stored in the inten- 
sity of the diffracted beams. As it is known, the diffracted 
intensities for crystals with defects (i.e., for nonperiodic 
structures) cannot be calculated with the formulas usually 
adopted for the structure factors V s of ideal crystals. How- 
ever, when an average periodic lattice can be defined, the 
total scattered intensity distribution in a kinematical (single 
scattering) approximation can be written as 15 
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FIG. 2. Experimental geometry used for strain measurements 
from LACBED patterns. Using the same illumination condition of 
Fig. 1, but tilting the cross-sectioned sample of about 3.5° and 
selecting the transmitted beam with the selected area aperture, the 
{333} Bragg contour become visible in the pattern and superim- 
posed to the image. The difference between the position of the 
Bragg contour in the silicon reference crystal and the defective 
layer allows to measure the strain along the direction parallel to the 
surface normal. 

where j. (s = 2 sin(0)/\) is the scattering vector, \(V g )\ 2 is 
equal to the square of the gth component of the Fourier 
transform of the electrostatic potential ( V(r)) of the average 
lattice, which produces the diffracted intensity ! g at recipro- 
cal lattice points. The quantity |(AV 5 (.y))| 2 is the square of 
the Fourier transform of the deviation from the average lat- 
tice (AV(r)), which gives rise to the nonperiodic back- 
ground intensity (diffuse scattering) l d , close to the recipro- 
cal lattice point g. 

For the description of the defective layer in hydrogen- 
implanted silicon it seems reasonable to neglect the effect of 
hydrogen atoms on the electron-scattering factor; then in the 
simplest approach, the damage can be thought of as formed 
by small clusters of silicon point defects. Taking into account 
the relaxation of the lattice surrounding each cluster, which 
affects a number of atoms larger than those of the cluster 
itself, the decrease of Bragg intensity can be expressed by 
the structure factor of the crystal, vf, multiplied by an at- 
tenuation factor, 



(V g )=Vftxp[-L G l 



(2) 



For clusters having a relaxation volume Af} cj , present in an 
atomic concentration C d , the L G factor is given by 9,10 



L c =v5C cl g M 



(|Afl cl |) 3 ' 2 



(3) 



where g is the modulus of the reciprocal lattice vector and 
n Si is the silicon atomic volume at the equilibrium density 
(2X 10" 2 nnrVatom). The quantity |Aft c ,| depends on the 
degree of correlation of point defects forming the cluster, 
being equal, in a first approximation, to the atomic volume 
for noncorrelated point defects and to the volume of the 



slipped region for dislocation loops, respectively. The simple 
relation between the attenuation factor L G and the structural 
features of the damage reported in Eq. (3) is obtained assum- 
ing a linear superposition of the displacement field produced 
by each cluster present in the defective layer 9,10 In the case 
of point defect clusters intermediate between noncorrelated 
point defects and dislocation loops, Eq. (3) is assumed as a 
reasonable approximation of the effect of lattice disorder on 
diffracted beam intensity. 

In electron-diffraction experiments the reduction of the 
diffracted intensity due to the L G factor cannot be predicted 
on the basis of a simple kinematical approximation because 
of the well-known many-beam (multiple-scattering) charac- 
ter of the electron-diffraction process. However, in a previ- 
ous paper 12 it has been found that when the sample is ori- 
ented along a high indices zone axis, the ratio R between the 
integrated intensity in the silicon substrate, if, and in dam- 
aged silicon, I g , is given by 



R = 



= exp[26L c ]. 



(4) 



In particular, for the (880) diffracted beam of the [44 1] zone 
axis the coefficient b was found to be equal to 1.26 for 225- 
330 nm thickness of the sample, which is useful for 
LACBED experiments at 100 keV. The constant b defines the 
correction factor, which is applied to the logarithm of the 
intensity ratio reported in Eq. (4) in order to determine L G in 
a "quasikinematicar approximation. Then the experimental 
method for disorder depth profile using DFLACBED pat- 
terns is simply based on the measurement of the profile R(z) 
of the ratio R reported in Eq. (4), taking the integrated inten- 
sity in the silicon substrate as a reference. Once L G (z) is 
determined, it is possible to deduce Af2 cl and C c , if the local 
strain is determined. 

B. Strain measurements 

Local strain analyses inside the defective layer have been 
performed following the method proposed for Si/SiGe 
heterostructures, 13 which supplies a simple relationship be- 
tween lattice strain and bending of Bragg contours observed 
in BFLACBED patterns when they cross a strained layer. 
Here this method is preferred to the high-order Laue zone 
(HOLZ) lines method used in heterostructures 14,16 because, 
due to the high level of static disorder present in the dam- 
aged layer, HOLZ lines become invisible, whereas the Bragg 
contours of medium index reflections, such as {333}, should 
be detectable even in disordered crystals. 

In the case of ion implanted materials, a biaxial lattice 
strain can be present. The perpendicular strain [s l = (a 1 
~ a si /a si)] evidences the change of the lattice parameter 
along the direction perpendicular to the surface (a L ) with 
respect to the undamaged substrate (a Sl ). The parallel strain 
O„ = (<2,|-0 Si )/a Si ] shows the variation of the in-plane lat- 
tice constant (a„) with respect to the silicon lattice along a 
direction parallel to the surface. The strain e„ is present when 
the implanted layer undergoes plastic deformation caused by 
the injection of stress- relieving defects. When both e L and e H 
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are not zero, they could be measured using two independent 
reflections from a couple of equations such as the following: 



2A<p (; '* /) 2q {hkl) /p 



(5) 



sin2<p sin2<p 

Here ti^ 0 is the strain induced variation of the angle <p 
between the (hkl) plane and the (001) interface planes, p is 
the .angular distance from the pole corresponding to the in- 
tersection of {hkl) and {hkl) contours, and q (hki) is the shift 
of the Bragg contour (see Fig. 2). 

When the <?„ component is negligible (an approximation 
frequently used in ion implanted silicon 17 ) Eq. (5) applied to 
the {333} Bragg contour of the [ I [0] cross section samples 
reduces to: 



2A<p 



,(333) 



sin2<p 



lq™lp 
sin 2 <p 



(6) 



Once e x is known, the stress in the plane of the interface can 
be determined using the following equation: 



Ev 



(J— 



(7) 



where E and v are the Young's modulus and the Poisson's 
ratio, respectively. Here we assume that the isotropic elastic 
constant of the damaged layer are, in a first approximation, 
equal to those of the silicon substrate [E= 130.2 GPa and v 
= 0.28 (Ref. 18)]. 

In the frame of the cluster model described in the preced- 
ing section, perpendicular strain can be related to both the 
cluster concentration and the relaxation volume through the. 
following relation: 9,10 ' 17 



1 4- v AH cJ 



"3(i-v) Ucl n Si ' 



(8) 



where the term ( 1 + v)l{ 1 - v) accounts for the uniaxial ge- 
ometry of the strain field. 

IV. RESULTS 

A. Electron microscopy 

In Fig. 3 the TEM images of the analyzed samples are 
reported. The as implanted sample [Fig. 3(a)] shows a defec- 
tive region centered at a depth of 210 nm and a width of 80 
nm wide. Due to the small dimensions (diameter of 3 nm) 
and to the high number density («* 10 12 cm" 2 ) of the defects 
it is difficult to characterize them. Rising the annealing tem- 
perature (300 °C, 2 h), platelets (diameter of 12 nm, number 
density of 5X 10 11 cm" 2 ). become clearly visible [Fig. 3(b)]. 
Lattice images of the defective layer (not reported here) al- 
low to measure the thickness («*l nm) and to deduce the 
zigzag morphology of platelets having {111} and {100} mean 
habit planes. The out-of-focus images of platelets show a 
clear' variation in contrast from underfocus to overfocus, 
typical of defects having an electron-optical density lower 
than that of the silicon matrix, thus suggesting hydrogen 
decoration. 2 The analysis of the defect distribution inside the 
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FIG. 3. Cross section TEM images of the analyzed 
samples, (a) Dark field (£ = 004) image of the as implanted 
sample, (b) Out-of-focus image of the 300 °C 2 h annealed sample. 
{111} and {100}: platelets are visible, (c) Out-of-focus image of 
the 500 °C 2 h annealed sample. Some "microcracks" with {111} 
and {100} facets are detected. 

defective layer shows that {111} platelets are mainly located 
in the deep tail of the band thus suggesting a negligible role 
of the implantation damage on their formation. This obser- 
vation is in agreement with experimental results obtained on 
plasma hydrogenated silicon samples, where, in absence of 
implantation damage, {111} platelets form. 19 On the contrary, 
{100} platelets are detected mainly in the central part of the 
defective layer where the radiation damage and the hydrogen 
profiles mainly intersect. Dimensions and concentration of 
platelets are in agreement with literature data on the study of 
the growth kinetics at 450 °C for 1-15 min in samples im- 
planted at room temperature at a dose of few 10 16 cm* 2 that 
have shown a conservative ripening (Ostwald ripening) of 
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the extended defects. 20 In Fig. 1(c) the TEM cross-section 
image of the 500 °C annealed sample is reported; extended 
defects have lateral dimensions and a number density similar 
to those observed in the 300 °C annealed samples but show 
{100} and {111} faceted cores. On the basis of this TEM 
evidence platelets observed after 300 °C annealing are distin- 
guished from defects observed in the 500 °C annealed 
sample, hereafter called "microcracks." 

B. LACBED static disorder and strain measurements 

In Fig. 4(a) is reported a (880) DFLACBED, aside of the 
as-implanted sample image in order to compare the intensity 
of the Bragg contour with the morphology of the sample. 
The diffracted intensity decreases when it reaches the defec- 
tive region and presents a minimum approximately at the 
depth where the image shows the center of the defective 
band (210 nm). 12 

The static disorder increase, which follows the annealing, 
is qualitatively observed by comparing the DFLACBED pat- 
tern of the sample annealed at 300 °C [Fig. 4(d)] to the as 
implanted one, as the (880) Bragg contour loses almost com- 
pletely its intensity in the damaged region of the annealed 
sample. An appreciable recovery of crystal quality is de- 
tected in the sample annealed at 500 °C [Fig. 4(f)]. The raise 
of the static disorder strength from the as implanted sample 
to the 300 °C annealed one can be better analyzed if the 
depth profiles of the integrated intensity ratio R(z) of Eq. (4) 
measured in the three analyzed samples are compared. It 
must be pointed out that in the defective layer of the an- 
nealed samples of the present study, due to the high level of 
static disorder, the diffracted beams loose a large part of their 
intensity and high diffuse scattering is observed near the 
Bragg position. For this region, we have tried to improve the 
process to give a more reliable method of disorder depth 
profiling proposed in Ref. 12, thereby extending it to the case 
of highly damaged samples. The procedure proposed here 
differs from the previous one in two aspects: first, the inte- 
grated intensity profile is directly obtained from the 
DFLACBED patterns and not from the envelope of single 
rocking curves. Second the background subtraction is per- 
formed using depth profiles of integrated diffuse intensity 
and not using a linear background subtraction on each rock- 
ing curve. In Fig. 5 the steps followed in order to obtain the 
depth profile of the integrated intensity are reported. 

(i) Directly from the DFLACBED [Fig. 5(a)] three inte- 
grated intensity profiles are obtained. The first one, PI, al- 
lows to obtain the raw integrated intensity profile of the 
Bragg contour [Fig. 5(b)]. The profiles P2 [Fig. 5(c)] and P3 
[Fig. 5(d)] are used to perform the background subtraction. 

(ii) Profiles P2 and P3 are averaged [Fig. 5(e)] in order 
to approximately evaluate the integrated background under, 
the Bragg peak. In principle, this point is the most critical 
one but seems reasonable because, as expected for point de- 
fects and clusters, 21 the peak of the diffuse scattering is 
asymmetric and not coincident with the Bragg one. 

(iii) The average profile obtained in step (ii) is subtracted 




FIG. 4. Cross section TEM images and (880) DFLACBED pat- 
terns of the as implanted sample (a), (b); 300 °C 2 h annealed 
sample (c), (d); and 500 °C 2 h annealed sample (e), (f). The images 
allow to calibrate the depth scale of the DFLACBED patterns. 

to P 1 in order to obtain the background subtracted integrated 
intensity depth profile / ? (z) [Fig. 5(f)]. 

In Fig. 6 are reported the R(z) plots obtained by applying 
Eq. (4) to the integrated intensity profiles after background 
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FIG. 5. Procedure adopted in order to obtain the depth profile of the integrated intensity (intensity), (a) DFLACBED pattern of the 300 °C 
2 h annealed sample. The rectangular masks correspond to the integration field of the rocking curve, near the Bragg peak (PI), and in the 
diffuse ^scattering region close to the Bragg peak (P2) and (/>3). (b) Raw diffracted intensity profile, PI. (c) Diffuse intensity profile po 
(d) Diffuse intensity profile, P3. (e) Average diffuse intensity profile, (f) Background subtracted intensity profile 
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FIG. 6. Depth profile of the ratio between the integrated inten- 
sity in the defective layer and in the silicon substrate as a function 
of the annealing temperature. 

subtraction and calibration of the depth scale. Now the quan- 
titative comparison between the crystal quality of three 
samples can be easily done. 

Strain measurements have been performed using the 
LACBED technique in the BFLACBED geometry (see Sec. 
DIB). In order to quantify the perpendicular strain e ± the 
angular shift q of the Bragg contour corresponding to the 
{333} crystallographic planes inside the defective layer is de- 
termined. Actually, due to the well-known relaxation of 
stresses in cross-sectioned samples, the angular shift q comes 
out from the unrelaxed part of the strain field. This system- 
atic error is of the order of 15% in (110) cross sections 16 and 
is taken into account in the evaluation of the error bar of the 
strain values. In Fig. 7(a) the experimental LACBED pattern 
of the 300 °C annealed sample that allows to measure the 
shift of the (333) reflection inside the damage region is 
shown. Due to the strong static disorder the Bragg contour is 
barely visible in the image. In Fig. 7(b) an edge-filtered im- 
age of the {333} Bragg contour of Fig. 7(a) is reported, show- 
ing slightly better the angular shift of the Bragg contour 
inside the damaged region. In Fig. 7(c) line scans across the 
{333} reflection are reported, near the middle of the damage 
layer, showing well-defined minima that have been used to 
measure the angular shift q from the silicon Bragg position 
(reported for comparison). The angular shift of the {333} 
reflection in the 500 °C annealed sample is indistinguishable 
from the as-implanted one and is omitted in the plot. In fable 
I the mean values of both the perpendicular strain and the 
related in-plane stress are reported showing the high stress 
(compressive, about 3 GPa) present in the defective layer of 
the 300 °C annealed sample. A more accurate analysis of the 
stress/strain profile was impossible because of the inclined 
position of the {333} Bragg contour with respect to the dam- 
age layer, which worsens the spatial resolution. 

V. DISCUSSION 

From a TEM point of view, platelets (clearly visible in the 
300 °C annealed sample) seem to be the relevant extended 
defects formed after low-temperature annealing. After an- 




q (radians) 

FIG. 7. (a) BFLACBED taken in the 300 °C annealed sample 
for the strain measurement. The {333} and {444} Bragg contour are 
vertical and inclined of about 54° with respect to the surface. The 
bending of the Bragg contour inside the defective layer is visible. 

(b) Edge filtered image of the {333} Bragg contour to make more 
evident the angular shift of the Bragg contour in the defective layer. 

(c) Intensity profiles across the {333} Bragg contour used for strain 
measurement in the analyzed samples. 

nealing at 7*=500°C these defects evolve towards "microc- 
racks," i.e., faceted intrinsic defects with a nonequilibrium 
shape when, as reported in the literature, 2 * 3,7 H starts to des- 
orb from the sample. Electron-diffraction measurements of 
strain and static disorder show a "reverse annealing." In par- 
ticular, the 300 °C annealed sample appears more stressed 
and disordered than the other samples. These evidences are 
in qualitative agreement with previously reported x-ray 
double crystal diffraction DCD (Ref. 6) and RBSC (Ref. 2) 
measurements performed on similar samples. In particular, 
strain and static disorder measured by DCD and LACBED 

TABLE I. Mean perpendicular strain and in-plane stress values 
obtained from LACBED measurements in the three analyzed 
samples 



As implanted 300 °C 500 °C 

Sjl (10~ 3 ) .5±1 12±2 5±1 

a- (GPa) - 1 ±0.1 -3 ±0.2 -I ±0.1 
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FIG. 8. Plot of the mean relaxation volume and of the mean 
concentration of clusters as a function of the annealing temperature. 
The values are deduced by combining strain and static disorder 
results. 

are, in an isotropic approximation, almost directly compa- 
rable as diffraction techniques are mainly sensitive to the 
long-range displacement fields produced by defects. The 
agreement is satisfactory in as-implanted and 500 °C an- 
nealed samples, where DCD measurements show the pres- 
ence of a perpendicular strain almost constant inside the de- 
fective layer [ranging from (5-7) X 10~ 3 ] and the highest 
value of static disorder in the annealed sample. Slightly dif- 
ferent is the situation in the 300 °C annealed sample, where 
BFLACBED measurements show a perpendicular strain 
value higher than 10"" 2 , which has no counterpart in DCD 
measurements probably because the peak of static disorder is 
so high that, for x-ray measurements, the presence of an 
amorphous buried layer at the depth of about 210 nm cannot 
be excluded. 6 On the contrary, from electron-diffraction 
and high-resolution TEM (not reported here) point of view, 
this sample appears still crystalline, but characterized 
by a high level of static disorder and a measurable strain. 
Strain and static disorder data allow to exploit the cluster 
model reported in Sec. Ill for the quantitative structural 
interpretation of electron-diffraction measurements in 
terms of cluster mean relaxation volume and cluster concen- 
tration [see Eqn. (3) and Eqn. (8)]. In Fig. (8) the results 
obtained by applying the cluster model to the data measured 
near the peak of the damage are summarized. According 
to strain measurements, the mean relaxation volume is posi- 
tive, but always higher than that expected in the 
case of interstitial point defects (Af2,= 1.8X 10" 2 nm 3 
= 0.9 H Si 22 ). Taking the data from the as-implanted sample 
as a reference for comparison among the samples, it is found 
that the mean relaxation volume increases by a factor 3 in the 
300 °C annealed sample and by an order of magnitude after 
annealing at 500 °C. On the contrary the cluster concentra- 
tion does not change significantly upon the first annealing 
step and only after annealing at 7=500°C it markedly de- 
creases. It can be further observed that the relaxation volume 
(cluster concentration) determined by DFLACBED is always 
lower (higher) than the corresponding quantities deduced 
from extended defects observed by TEM. 2 These evidences 
support the idea that the extended defects represent only a 
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FIG. 9. Comparison between the concentration of displaced at- 
oms expressed in numbers of "interstitial equivalent point defect" 
as determined by LACBED and RBS ch. 

small fraction of the total damage distribution, whereas 
hydrogen-related point defect clusters grow upon low- 
temperature annealing. In addition, the increase of the relax- 
ation volume at constant cluster concentration observed after 
annealing at 300 °C suggests a nonconservative growth of 
the damage during this .thermal treatment. Then the experi- 
mental evidences reported in Fig. 8 could be the structural 
fingerprint of the progressive segregation of H 2 molecules in 
hydrogen decorated vacancies, which, following a recently 
proposed model, 8 sustain the growth of second-neighbor 
VH 4 's with the emission of interstitial defects, thus generat- 
ing very complex clusters of damage. If this attribution is 
correct, these clusters could have a relaxation volume of 
about 0.1 nm 3 . The additional relaxation volume increase 
and the significant cluster concentration decrease observed 
after 500 °C annealing agree with the model of conservative 
growth of, damage observed for extended defects in samples 
annealed at 450 °C. 20 

In order to proceed in the comparison among electron- 
diffraction measurements and results obtained with different 
techniques on the same samples, it has been tried to compare 
the electron-diffraction data, obtained from the cluster 
model, to RBSC results. Differently from the case of diffrac- 
tion (electron or x-ray) techniques, this comparison could be 
performed if electron-diffraction and ion techniques both 
convert "raw" data (strain and static disorder for electron 
diffraction and dechanneling yield for RBSC) in comparable 
quantities using models of damage. A quite common inter- 
pretation of RBSC data is to convert dechanneling yield 
depth profiles in "interstitial equivalent" displaced atoms 
profiles (see, e.g., Ref. 2). The same parameter can be ob- 
tained from electron-diffraction data, by simply multiplying 
the cluster relaxation volume, the cluster concentration, and 
the silicon atomic density (5 X 10 22 atoms/cm 3 ) and normal- 
izing to the single interstitial relaxation volume. Here the 
additional assumption with respect to the cluster model re- 
ported in Sec. Ill A is that not only the displacement field 
produced by each cluster, but also the displacement field pro- 
duced by each atom in the cluster superimpose linearly. In 
Fig. 9 are reported the defect concentrations measured by 
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DFLACBED near the damage peak as a function of the an- 
nealing temperature. The trend is found in fair agreement 
with the RBSC data reported in the literature" thus support- 
ing the "nonconservative" growth of damage in the low- 
temperature annealed sample. The agreement in the absolute 
values of the defect concentration is not satisfactory, particu- 
larly in the 500 °C sample where the dimension of defects 
make the hypothesis of linear superposition of the displace- 
ment field inside each cluster not applicable. 

VI. CONCLUSIONS 

In this paper the TEM/LACBED techniques have been 
applied to characterize the defective layer in hydrogen - 
implanted single-crystal silicon. It has been shown that the 
LACBED technique is able to measure strain and static dis- 
order in heavily damaged crystals, allowing to determine the 
mean relaxation volume and the mean concentration of clus- 
ters present in the defective layer. It is found that for samples 
annealed at 7=300°C for 2 h platelets (about 12 nm long 
and 1 nm thick) form, but smaller clusters of point defects 



exist in larger concentration, as the mean relaxation volume 
of defects is of the order of 0.1 nm 3 . In this sample the 
cluster mean relaxation volume is higher than that measured 
in the as-implanted sample, whereas the mean concentration 
of defects does not vary appreciably. An additional increase 
of the relaxation volume and a significant decrease of the 
number density of clusters is found after annealing at 500 °C 
for 2 h. These evidences indicate that only after this thermal 
treatment, responsible in the mean time for the hydrogen 
effusion from the surface and for the formation of "microc- 
racks," defects seem to start the surface reconstruction and 
the conservative ripening, which, in samples implanted at a 
higher hydrogen dose, is the base of the "smart-cut" techno- 
logical process. 
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